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Abstract A statistic based on increment ratios (IR) and related to zero crossings of increment sequence is defined 
and studied for measuring tlie roughness of random paths. The main advantages of this statistic are robustness to 
smooth additive and multiplicative trends and applicability to infinite variance processes. The existence of the IR 
statistic limit (called the IR-roughness below) is closely related to the existence of a tangent process. Three particular 
cases where the IR-roughness exists and is explicitly computed are considered. Firstly, for a diffusion process with 
smooth diffusion and drift coefficients, the IR-roughness coincides with the IR-roughness of a Brownian motion and 
its convergence rate is obtained. Secondly, the case of rough Gaussian processes is studied in detail under general 
assumptions which do not require stationarity conditions. Thirdly, the IR-roughness of a Levy process with a— stable 
tangent process is established and can be used to estimate the fractional parameter a £ (0, 2) following a central 
limit theorem. 
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oo ! 1. Introduction and the main results 
O 

It is well-known that random functions are typically "rough" (non-difFerentiable), which raise the question 
, of determining and measuring roughness. Probably, the most studied roughness measures are the Hausdorff 

H ' dimension and the p— variation index. There exists a considerable literature on statistical estimation of 

these and related quantities from a discrete grid. Hence, different estimators of the Hausdorff dimension 
have been studied, as the box-counting estimator (see Hall and Wood, 1993 for stationary Gaussian processes 
or Levy-Vehel and Peltier, 1994, for Gaussian processes with stationary increments). To our knowledge, the 
_ff-variation estimator, where is a measurable function, was first proposed by Guyon and Leon (1989) 
for stationary Gaussian processes where central and non-central limit theorems are established following 
the Hermite rank of H and the asymptotic local properties of the variogram and its second derivative. 
Further studies provided a continuation of this seminal paper in different ways. Istas and Lang (1997) 
studied generalized quadratic variations of Gaussian processes with stationary increments. CoeurjoUy (2001 
and 2005) studied £P-variations of fractional Brownian motion and £^-variations of multifractional Brownian 
motion. CoeurjoUy (2007) discussed L-variations based on linear combinations of empirical quantiles for 
Gaussian locally self-similar processes. An estimator counting the number level crossings was investigated 
by Feuerverger et al. (1994) for stationary Gaussian processes. 



*This author thanks the Steklov Institute of St Petersburg (Russia) for a nice and fruitful invitation 
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In the present paper we introduce a new characteristic of roughness, defined as a sum of ratios of consecutive 
increments. For a real-valued function / = {f{t),t G [0, 1]), define recursively 

Aj'"/ /(^) -/(-)- 
n n 

Af-f A]^"Af 1'"/ = j^i-^r^'ifjfC-^)^ (1-1) 

i=o ^^"^ ^ 

so that A^'"/ denotes the p— order increment of / at ^, p = 1, 2, • • • , j — 0,1, ■ ■ ■ ,n — p. Let 

1 I AP^"f -4- A^'" fl 



with the convention ^ := 1. In particular, 



-2 



" --itoi/(^)-/(^)i+i/m-/(^)r 

Note the ratio on the right-hand side of (ll.2p is either 1 or less than 1 depending on whether the consecutive 
increments A^.'"/ and A^^j^/ have same signs or different signs; moreover, in the latter case, this ratio 
generally is small whenever the increments are similar in magnitude ("cancel each other"). Clearly, < 
< 1 for any f,n,p. Thus if limi?P'"(/) exists when n — ^ oo, the quantity can be used to 

estimate this limit which represents the "mean roughness of /" also called the p—th order IR-roughness of 
f below. We show below that these definitions can be extended to sample paths of very general random 
processes, e.g. stationary processes, processes with stationary and nonstationary increments, and even L'- 
processes with q < 1. 

Let us describe the main results of this paper. Section[2]derives some general results on asymptotic behavior 
of this estimator. Proposition 2.1 says that, for a sufficiently smooth function /, the limit lim„^oo R^'"{f) = 
1. In the most of the paper, f = X is a random process. Following Dobrushin (1980), we say that 
X = {Xt, t e M) has a small scale limit at point io G K if there exist a normalization A'^*°'>{S) oo 

when (5 — >■ and a random process 

y(to) = (!;(*«), r > 0) such that 
A^'"\6)iXt,+rS-Xt„) (1.4) 

where -L^:^' stands for weak convergence of finite dimensional distributions. A related definition is given 
in Falconer (2002, 2003) who called the limit process a tangent process (at to). See also Benassi et 

al. (1997). In many cases, the normalization A^*°\S) = 6-^^*°\ where < i?(io) < 1 and the limit tangent 
process is self-similar with index -ff(to) (Falconer, 2003 or Dobrushin, 1980). Proposition 2.2 states 

that if X satisfies a similar condition to (|1.4I) . then the statistic W'^{X) converges to the integral 



i?P'"(X) A / E 

n-i-oc Jo 



|Agy(*) + APyW| 
|Aj;r(*)| + |APr(*) 



At, (1.5) 



in probability, where A^F^*^ = A^'^F*^*^ — X^f^ol"!)^ J = 0; 1 is the corresponding increment 

of the tangent process F^*-' at t G [0, 1). In the particular case when X has stationary increments, relation 
(11.51) becomes 



RP'"{X) E 



(1.6) 



|Agy| + |Afy| 

Section [3] discusses the convergence in (jl.Sp for diffusion processes X admitting a stochastic differential 
dX = atdB{t) + btdt, where i? is a standard Brownian motion and (at), (6t) are random (adapted) functions. 
It is clear that under general regularity conditions on the diffusion and drift coefficients (at), (bt), the process 
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X admits the same local Holder exponent as B at each point t^ G (0, 1) and therefore the IR-roughness 
of X in p. 51) should not depend on these coefficients and should coincide with the corresponding limit for 
X = B. This is indeed the case since the tangent process of X at t is easily seen to be 1^'*^ = otB and the 
multiplicative factor at cancels in the numerator and the denominator of the fraction inside the expectation 
in (|1.5|) . See Proposition 3.1 for details, where the convergence rate 0{n^^^) (a.s.) in (jl.5|) with explicit limit 
values Ap(l/2) is established for diffusions X and p = 1,2. 

Considerable attention is given to the asymptotic behavior of the statistic RP''^{X) for "fractal" Gaussian 
processes (see Section 14]). In such a frame, fractional Brownian motion (fBm in the sequel) is a typical 
example. Indeed, if X is a fBm with parameter H G (0, 1), then X is also its self tangent process for any 
t € [0, 1] and (see SectionE]): 



RP'^^iX) 



Ap{H) {p=l,2), 



(1.7) 



where 



Ap(ff) 



HPpiH)), 



A(r) := — arccosi— r)H — \ log, ^ 

^ ' TT ^ ' 7rVl-?' Vl+7' 



PviB) 



coTT{Af,BH,APBH) 



(1.8) 
(1.9) 
(1.10) 



and where A]Bh = BrU + 1) - BhU), A^Br = BhU + 2) - 2BhU + 1) - BhU) (j € Z) are respective 
increments of fBm. Moreover, 



V^(i?f'"(X)-Ap(i/)) A AA(0,Ep(i7)) if 

n— >-oo 



p=l, 0<H<3/4, 
p = 2, 0<H <l, 



(1.11) 



where — > stands for weak convergence of probability distributions. The asymptotic variances T,p{H) in 
(|l.lip are given by 



S,(i/) :=^cov ^ 







\^oBh\ - 





(1.12) 



The graphs of Ap{H) and y/'Sp{H) {p — 1, 2) are given in Figures [T] and [2] below. 





Figure 1. The graphs of Ai{H) (left) and A2(H) (right). 



The difference in the range of the parameter H for p = 1 and p = 2 in the central limit theorem in (II. lip 
are due to the fact that the second order increment process (A|_B^f,j e Z) is a short memory stationary 
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Figure 2. The graphs of ^T,p{H), p = I (with a pole at 3/4) and p = 2 (with a pole at 7/4) (from Stoncelis and Vaiciulis, 
2008, with kind permission of the authors) 

Gaussian process for any H e (0, 1), in contrast to the first order increment process (Aji?jif, j G Z) which 
has long memory for H > 3/4. 

Generahzations of (|1.7p and (|l.lip to Gaussian processes having nonstationary increments are proposed 
in Section m Roughly speaking, RP''^{X), p = 1,2 converge a.s. and satisfy a central limit theorem, 
provided for any t e [0, 1] the process X admits a fBm with parameter H(t) as a tangent process (more 
precise assumptions (A.l), (A.l)' and (A.2)p are provided in Section HJ. In such frames, the limits in 
()1.7p are Ap{H(t))dt instead of Ap{H) and the asymptotic variances in also change. The case 

of Gaussian processes with stationary increments is discussed in detail and the results are used to define 
a -^71" consistent estimator of H, under semiparametric assumptions on the asymptotic behavior of the 
variogram or the spectral density. Bardet and Surgailis (2010) study a punctual estimator of H{to) obtained 
from a localization around to £ (0, 1) of the statistic i?^^"(X). 

The main advantages of estimators of the type (jl.2p involving a scaling invariant function of increments 
seem to be the following. Firstly, the estimator i?P'"(X) essentially depends on local regularity of the process 
X and not on possible "multiplicative and additive factors" such as diffusion and drift coefficients in Section 
[3] or smoothly multiplicative and additive trended Gaussian processes, see Proposition l4.1l of Sectional This 
property is important when dealing with financial data involving heteroscedasticity and volatility clustering. 
Such a robustness property (also satisfied by the estimators based on generalized quadratic variations of 
wavelet coefficients) represents a clear advantage versus classical parametric Whittle or semi-parametric log- 
periodogram estimators. Secondly, the estimators in (|1.2p are bounded functionals and have finite moments 
of any order. Section [5] discusses jump Levy processes, with the Levy measure regularly varying of fractional 
index a € (0,2) at the origin. Using a modification of (|1.2p . we define a -^/n— consistent estimator of a, 
together with a central limit theorem, in a very general semiparametric frame. This result is new and 
interesting because there exist very few papers providing consistent estimators of a (to our knowledge, the 
only comparable results have been established in (Belomestny, 2010) and (Ait Sahalia and Jacod, 2009) in 
a financial and somewhat different context). Finally, in the Gaussian case, using the approximation formula 
provided in Remark 14. 3[ an estimator of H based on i?^'"(X) can be extremely simply computed: 

^ ^ / \ V— ^ L\^ fc + 2 — 2X k + \ + X k_ -\- X k + -^ — 2X k + 2 -\- X k + 1 

jj(2) ^ / I 2 2 2^ 2 2 2_L _ 5174 

0.1468 Vn - 2 ^ - 2Xii+i + I + - 2Xfc+2 + 1 

In the R language, if X is the vector {Xi, , X2 , • • • , Xx) , 

H';'^^ ~ (mean(abs(diff(diff(X[-l])) + diff(diff(X[-length(X)]))) 

/(abs(dif f (dif f (X[-l]))) + abs(dif f (dif f (X[-length(X)]))))) - 0.5174)/0.1468. 



Therefore its computation is very fast and does not require any tuning parameters such as the scales for 
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estimators based on quadratic variations or wavelet coefficients. The convergence rate of our estimator is ^Jn 
as for the parametric Whittle or the generalized quadratic variation estimators and hence it is more accurate 
than most of other well-known semi-parametric estimators (log-periodogram, local Whittle or wavelet based 
estimators). 

Estimators of the form (|1.2|) can also be applied to discrete time (sequences) instead of continuous time 
processes (functions). For instance Surgailis et at. (2008) extended the statistic i?^'"(X) to discrete time 
processes and used it to test for I{d) behavior (—1/2 < d < 5/4) of observed time series. Vaiciulis (2009) 
considered estimation of the tail index of i.i.d. observations using an increment ratio statistic. 

Remark 1.1 The referee noted that the IR-roughness might be connected to the level crossing index (see 
Feuerverger et at, 1994). To our surprise, such a connection indeed exists as explained below. Let Yn{t), t € 
[0, 1 — ;^] be the linear interpolation of the "differenced" sequence A^-'"X — X {■^^) — X {^) , j ^ 0,1, ■ ■ ■ , n—1: 



(^-t)A]-"X + (t-^)A]-X 



Yn{t) = n 

j = 0,1, ■ ■ ■ ,n — 2. Then, using Figure [3] as a proof, 



n n 



i?i."(X) = -^Y\^(=^s\te[-,^-^):Y^{t)>o]-mcas\te[^,^-^):Y^{t)<o] (1.13) 
n — 1 ^-^ linn J Inn J 



n-2 





Figure 3. The proof of l|1.13|l : follows by 



\Y„(i-)+Y„(i±^)\ 



^)l 



n|(7i -(72| 



Let 'ip{xi,X2) '■— \xi -|-X2|/(|xi| -I- |a;2|), ^'0(2^1, 2:2) l(a;ia;2 > 0). Clearly, the two quantities 1 — 
4'{Yn{^),Yn{^i^)) and 1 — 'ipQ(Yn{^),Yn{^^)) both are strictly positive if and only if y„ crosses the zero 
level in the interval '^^) but the former quantity measures not only the fact but also the "depth" of the 
crossing so that 1 — ■!/)(y„(^), Yn{^^^)) attains its maximal value 1 in the case of a "perfect" crossing in the 
middle of the interval ^^-^) (see Figure 

It seems that similar asymptotic results can be obtained for i?Q'"(X) :— ^JJ^p"^ i/)o(A^'"X, A^;|\X) 
measuring the number of zero crossings of the increment sequence A^'"X, k — 0, 1, • • • ,n — p and other 
similar statistics obtained by replacing the functions tp or ^0 by other scaling invariant functions. Let us 
note that i?Q'"(X) is related to the zero-crossings' counting statistic studied in Ho and Sun (1987) for sta- 
tionary Gaussian time series. Also note that the Hcrmite rank of ■(/'o is 2 and that the corresponding limit 
function Ao(r) = arccos(— r) is strictly increasing on the interval (—1, 1) similarly as the function A(r) in 
(|1.9p . On the other hand, while the statistic R^'^{X) is certainly of interest, the statistic RP'"'{X) seems 
preferable to it for the reasons explained above. In particularly, in the case of symmetric Levy processes X 
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with independent increments studied in Section [Sj the latter statistic leads to an estimator of the fractional 
index while the former statistic can be easily shown to converge to 1/2. 

The paper is organized as follows. Section [5] discusses some general (consistency) properties of the esti- 
mators RP "{X). Section [3] deals with the case when X is a diffusion. The case of Gaussian processes X is 
considered in Section |4] while the case of Levy processes is studied in Section [5l Section |6] contains proofs 
and other derivations. 

Below, we write C for generic constants which may change from line to line. 



2. Some asymptotic results 

The definition of in (|1.2p can be extended to more general increments (the so-called generalized 

variations). Consider a filter a :— {ao, ■ ■ ■ , aq) E M'+^ such that there exists p E < q satisfying 

^tae^O for i = 0,---,p-l and ^Fa^T^O. (2.1) 

The class of such filters will be denoted A{p, q). For n G N* {1, 2, • • •} and a function / ~ {f{t), t G [0, 1]), 
define the generalized variations of / by 

A;-'7 J=0,l,---,"-g. (2.2) 

A particular case of (|2.2p corresponding to q = p > 1, ae = (— 1)p~^(^) is the p~order increment A^'"/ in 
prri) . For a filter a e A{p,q), let 

n — q—1 \ \a.n j* i \a.n r\ 

■ ^-1 h iAr /I + iA^v"i/r 

It is easy to prove that R^'"'{f) — > 1 if / is continuously differentiable on [0, 1] and the derivative /' 

does not vanish on [0, 1] except maybe for a finite number of points. Moreover, it is obvious that i?^'"(/) = 1 
if / is monotone on [0, 1]: the IR-roughness of a monotone function is the same as of a smooth function, 
which is not surprising since a similar fact holds for other measures of roughness such as the p— variation 
index or the Hausdorff dimension. 

We conjecture that i?P'"(/) ^ 1 and R°-"{f) 1 for any g > p > 1, a S A{p,q) and / : [0, 1] ^ M which 
is (p — 1) times differentiable and the derivative f^P^^^ has bounded variation on [0, 1] with the support 
supp(/(P~i)) = [0,1]. However, we can prove a weaker result. 

Proposition 2.1 Let f be [p — l)-times continuously differentiable {p > 1) with f^P^^'^ being absolutely 
continuous on [0, 1] having the Radom-Nikodym derivative g = {f^^^^'^) ■ Assume that g a.e. in [0, 1]. 
Then RP '^lf) — > I and R°-'\f ) — > 1 for any a e A{p,q), q > p. 

Proof. We restrict the proof to the case p = 2 since the general case is analogous. Using summation by 
parts, we can rewrite A°'"/ as 

A;^"/ = Efe.A?-/, (2.4) 

i=0 

where bi := I]fe=o J2e=o = 0; 1; ' ••7 9, ^<?-i = bq = and b I]Lo 5 Z^Li ^^'^^ 7^ ^ in view of the 
assumption a G A{2,q). 
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Assume n is large enough and for a given t £ (0, 1), let A:„(t) £ {0, • • • , n — 2} be chosen so that t £ 
[kn{'t)/n, {kn{t) + l)/n), therefore kn{t) = [nt] — 1. We claim that for a.e. t £ (0, 1) 



lim n Al'" J = bg{t), lim n A^''' / = bg{t). 
Using the fact that the function {X1.X2) ^ \ xi\^xi \ continuous on M^\{(0,0)}, we obtain 

. \hg{t) + hg{t)\ 



(2.5) 







+ 





n^oo \bg{t)\ + \bg{t)\ 



(2.6) 



for a.e. t £ (0, 1), where we used the fact that bg{t) 7^ a.e. Since for n > q, R"''^{f) can be written as 
R°''^{f) ~ Jq h°-'^{t)dt, relation R°-'^{f) — > 1 follows by the dominated convergence theorem and the 

fact that < /i'''"(t) < 1. 

Relations ()2.5p can be proved using the Lebesgue-Vitali theorem (see Shilov and Gurevich, 1967, Ch. 4, 
§10, Theorem 1), as follows. Consider the signed measure /x on Borel subsets of [0, 1/2]^ given by 



fi{A) = / g(a;i + a;2)dxida;2- 
J A 



Note Al"f = /i((A:/2n, (fc + 2)/2n] x (fc/2n, (fc + 2)/2n]), (fc = 0, 2). Since rectangles + 

/i] X [a::2,2^2 + h], {0 < Xi < Xi + h < 1/2, i — 1,2) form a Vitali system on [0,1/2]^, the above mentioned 
Lebesgue-Vitali theorem implies that 



) fc„(ti)+2 



2n 



fc„(t2) fc„(t2) + 2 
2n ' 2n 



(2.7) 



a.e. in [0, 1]^. Taking into account the form of the measure ^ and the limiting function in (12.71) . it follows the 
convergence r^^A^'"/J^/ = 0„(i,t) — > g(t) a.e. on [0, 1]. Next, for any fixed i — 0,1, ■ ■ the sequence of 

rectangles (^^^, x (^4^^, ^^M±!±H], „ = 1, 2, • • • is regularly contracting to (ti, ^2) e (0, 1)^ 

in the sense of (Shilov and Gurevich, 1967, Ch. 4, §10). Hence, using the lemma on p. 214 of the above 
monograph, it follows that nV((^^, ^ {t^iMll^ fc„(^.H.+2 ]) _^ ^(u±*^) a.e. in[0,lp, 

implying 



2 A 2,n 

n A 



k (t)+if — ^ SC*) a.e. on [0,1], for any i = 0,l,' 



Together with (j2.4p . this proves (j2.5p and the proposition. 



□ 



Let us turn now to the case when f{t) — Xt, t £ [0, 1] is a random process. Now and in all the sequel, 
RP'^{X), R°-'^{X) are denoted RP'^, respectively. Below we formulate a general condition for the 

convergence of i?^'" and i?"'" to a deterministic limit. 



Assumption (A): For a.e. pairs (^1,^2) G (0, l)^,ii 7^ ^2, for i = 1,2 there exist: 

(i) normalizations A'*'^(^) — 00 (^ — 0), 

(ii) (mutually) independent random processes = (y(**)(T),T £ [0, 1]), 
such that for S —?' 0, si ti, S2 t2 

[A('^\6){X^,+Sr~X^J,A^'^\5){Xs,+Sr-Xs,)) (y(*^)(r),y(*^)(r)). (2.8) 

Remark 2.1 Relation ()2.8p implies the existence of a joint small scale limit (y(*i) y(*2)') g^_g_ pg^jj. 
(^1,^2) G (0, 1), with independent components Note Assumption (A) and Proposition 12.21 below 

are very general, in the sense that they do not assume any particular structure or distribution of X. 
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Proposition 2.2 Let a = (ag, - ■ ■ ,aq) G A{p,q), 1 < p < q be a filter and let X satisfy Assumption 
(A). Assume m addition that P(|A^y(*)| > 0) = 1, j = 0,1 for a.e. t e (0,1), where A'^Z = A°^^Z = 
Z^Lo + Then 



E R" 



E 







|Agr(*)| + |A°y(*)| 



dt 



Proof. The statement follows from 



n-i-oo Jo 



|AgF(*)| + |A!fr(*) 



and E(i?°^"-Ei?'''")^ 



(2.9) 



(2.10) 



Write Ei?-^'" = /„ E/i;;^"(t)dt, where (c.f. ([M 



A(*)(l/n) ( 






) 


AW(l/n) ( 




+ 




) 


|A[5y(*) + 


Ajr(*)| 




/iy(i), 




|Agy(*)| + 


|A?r(*)| 







for a.e. t £ (0, 1), according to Assumption (A) and the continuous mapping theorem. Whence and from 
< ft,^" < 1 using the Lebesgue dominated convergence theorem the first relation in ()2.10p follows. Moreover, 



E(i?'^>" - Ei?"'")2 



\n — Q/ 



1 fi 



E[/i^"(i)/i;^"(t')]dMi' - (Ei?"^'")', 



q' JO JO 

and with the same arguments as previously and the independence of F*-*-* and F*^* •* when t ^ t', 



E 



|Agr(t)| + |A^r(*)| |Agy(t')| + |A°r(*' 



and therefore J^' E[hT{t)hT{n]dtdt' [j^ Eh'^^-{t)dt 

relation in ()2.10p and the proposition. 



= E/i^(t) •E/i5.(f) 

0, thereby proving the second 

□ 



The following easy but interesting corollary can also be added to this result. It proves that smooth multi- 
plicative or additive trends do no change the L^-asymptotic behavior of Let €^[0, 1] denote the class 
of all p— times continuously differentiable functions on [0, 1]. 

Corollary 2.1 Let a e A{p,q) and X satisfy the conditions of Provosition \2.2\ with A'^^\5) = 0{5~^) 
(5 0) for each t G [0, 1]. Assume that a G C^[0, 1], (3 G ^[0, 1], inffgfoaj a{t) > and supjg[o,i] 1^(^)1 < °° 
a.s. Define Z such that Zt = a{t) Xt + I3{t), t G [0, 1]. Then (2^ holds with RP''' = RP'''{X) replaced by 
RP'"'{Z). 

Proof. We consider here p >2 but the case p = 1 can be easily obtained. With a(^^) = a(^) + ■^a'(-|;) (l + 
o(l)) and Al^'"/3 = 0(;^), for a.e. t G (0, 1) and k = kn{t) as defined in the proof of Proposition 12.11 we 
deduce that 

Arz = «(^)Arx + ia'(^)A^'"X + o(i) sup \Xit)\+0{^), 
n n n n te(o,i) 

with a' = {jaj)o<j<q G A{p ~ 1,?)- Therefore, 

^'*n^)A:-)Z = a(M!))A(0(i)A-)X + io,(A(*)(i)A:>X)+o,(i^ 

a(i)Agr(*). 



V 
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In a similar way, for a.e. pairs {t,t') e (0,1)^, t ^ t' , we can verify the joint convergence in distribu- 
tion of r.v.'s (i)A-"(,)Z, (i)A-"(,)^,Z, A(0(i)A- A(*')(i)A:- to the Hmiting r.v.'s 
a(t)Agy(*),a(t)A'j'y(*),a(i')Agr(*'),Q(t')Ajr(*'). Now, the statement of the corollary follows by the argu- 
ment at the end of the proof of Proposition 12.21 □ 

Remark 2.2 By definition, the statistics i?^'" and i?"'" for a G A{p, q), 1 < p < q are invariant with respect 
to additive polynomial trends of order less than p; in particular, i?^'" is insensitive to a quadratic trend 
while i?^'" does not have this property. On the other hand, Corollary 12.11 (see also Proposition 14. ip says 
that under weak additional conditions on X , any sufficiently smooth additive or multiplicative trends do not 
affect the limit of as soon as p > 1. In the important special case when the limit process F*^*^ — Bh in 
Assumption (A) and (j2.9p is a fractional Brownian motion with parameter H G (0, 1) independent of t, the 
statistic RP'^ converges in mean square to the expectation E | jftf'j^i ^ |a | ~ ^iPpi^))' ^-f- (|1.8p - (|1.10p . 
Numerical computations show that the correlation coefficient Pp{H) is a monotone function of H for any 
p > I and tends to constant value —1 on the interval (0, 1) as p increases. Therefore, for larger values of 
p, the range of X{pp{H)) is rather small and i?^'" seems less capable to estimate H. A final reason for our 
concentrating on the "lower-order" statistics i?^'"', p — 1,2 in the rest of the paper is the fact that i?^'" 
satisfies the central limit theorem in (|l.lip on the whole interval H G (0, 1). 

3. Diffusions 

Let ^ ^ 

Xt = Xa+ [ a,dB(s) + / 6,ds, t e [0, 1] (3.1) 
Jo Jo 

be a diffusion (or Ito's) process on M. In p.ip . we assume the existence of a right-continuous filtration T = 
{J-f, t G [0, 1]), a standard Brownian motion B adapted to moreover, Os, bs, s G [0, 1] are adapted random 
functions satisfying \bs\ds < oo, a^ds < oo a.s., and Xq is a J^o^ measurable r.v. Write Et[-] — E[-\J^t] 
for the conditional expectation. Let Ai(l/2) = A(/9i(l/2)) ~ 0.7206 and A2(l/2) = A(p2(l/2)) ~ 0.5881. 
The proof of the following Lemma |3. II is given in Annexe. 

Lemma 3.1 Let ip{xi^X2) ■= \xi + X2I/ (|a;i| -I- |a;2|) (a^i, X2 G M), and let Zi, i ~ 1,2 be independent Af{0, 1) 
r.v. 's. Then for any random variables £,1,^2, 

\E^{Zi+£i,Z2+£2)~^HZi,Z2)\ < 20 max (EC,')'/'. (3.2) 

Theorem 3.1 Assume the following conditions: there exist random variables Ki,K2 such that < Ki < 00 
a.s., and such that, for any sufficiently small h > and any 0<t<t + h< 1, the following inequalities 
hold, a.s.: 

\at\ > Ki, ^tbl+u < K2 and Et(at+h. - at)' < K2h. (3.3) 

Then 

i?'''"-Ap(l/2) = 0(n-i/3) a.s. (p=l,2). (3.4) 

Proof. We restrict the proof to the case p — 1 since the case p = 2 is analogous. For notational simplicity, 
assume that n is odd. Define 

\Al-"X + Al'^^X\ 

Vn{k) := i„ , ■!„ Vn{k) := ^k/niVnik)], fijk) := ijnik) - ?7„(fc) (3.5) 
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and correspondingly write = K + K, K ■= (n-l)-' Yl=l v'n{k), Ki {n~l)-' Eto"^' Vnim, 
K2 - 1)^' Eiro^^^' Vn{2k + 1). As {Vn{2k),T^2k+2)/n, fc = 0, . . . , (n - 2)/2) is a martingale difference 

sequence, so by Burkholder's inequality, 

/(«-2)/2 \^ 

and therefore 

00 00 
Y,P{\Ki\ > n-'/') < Cj2n^/'n-^ < 00, 

n— 1 n— 1 

implying i?^]^ = 0{n^^^^) a.s. A similar fact holds for -R^2- Thus, it remains to prove 

<-Ai(l/2) = 0(n-i/3) a.s. 



(3.6) 



Observe 



where 



7y;(fc)-Ai(l/2) - Ek/r. 



I^i(fc) + 6(fc) + ^2(fc)+6(fc)l 

i^i(fc)+a(fc)i + i^2(fc) + 6(fc)i 



E 



\Zi{k) + Z2{k)\ 



\Z^ik)\ + \Z2ik)\ 



Zi{k) 
Uk) 



r,l/2 



,1/2 



.V2Ai'"i?, Z2(fc) nV^A^- i?, 

»(fc+l)/n 



fc/n ^^k/n 
(k+2)/n 



(^_l)dS(s)+ni/2 / ^ds. 



{k+2)/n ^ 



/ (^_i)dS(s)+ni/M 

J(fc+l)/n ^a/c/n ^ ./(fc+lj/Ti 



-ds. 



'(fc+l)/n ^^fc/n ' -/(fc+lj/ri ^fc/n 

According to Lemma [3?T] above. Wn{k) — Ai(l/2)| < 36maxi=i,2 {F,k/n£,iik)) and therefore 

|i?;-Ai(l/2)| < 36max{(E,/„ef(fc))i/3 :i=l,2,fc = 0,l,...,7^-l}. 

Whence, (|3.6p follows from the following fact: there exists a r.v. K < 00, independent of n and such that 
for any 71 > 1, fc = 0, . . . , n — 1, z — 1, 2 



Ek/n&{k) < Kn-\ a.s. 



(3.7) 



Indeed, using (|3.3 

Efc/„^f(fc) = n 



{k+l)/n 

2n.^ _ . / E 



k/ni 1) 



1 ds + nEfc 



kjn ^O'k/ 
(fe+l)/ri 



k/n 



(k+l)/n ^ 



-ds 



k/n <^k/r: 



< nK 



-2 



'k/n 

'-2„-l 



Efe/„ (^fls - ak/n^ ds + 



{k+l)/n 



k/r. 



Efc/„6^ds 



and the bound (|3.7p for i = 2 follows similarly. This proves p. 71) and Theorem 13. 11 too. 



□ 



Let us present some examples of Ito's processes X satisfying conditions p.3p . 
Example 3.1 Let (X(, t € [0, 1]) be a Markov process satisfying a stochastic equation 

Xt = xo+ I a{Xs)AB{s) + / b{Xs)As, 

JQ Jo 



(3.8) 
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where G M is nonrandom, a{x), b{x), a; € R are real measurable functions and i? is a standard Brownian 
motion. Let Ft (j{B{s),s < t},0 < t < 1 be the natural filtration. Assume that 

\a{x)-a{y)\<K\x^yl \h[x) ~ b{v)\ < K\x - y\ (^,yeM) (3.9) 

for some constant K < oo. Then equation (j3.8p admits a unique adapted solution; see e.g. Gikhman and 
Skorohod (1969). Let at = a{Xt),bt — b{Xt). Assume in addition that \a{x)\ > Ki, {x e M) for some 
nonrandom constant Ki > 0. Then the first inequality in p.3|) is trivially satisfied; moreover, the second 
and third relations in p.3|) are also satisfied, with K2 = C(l + supQ<t<]^ Xf) < 00 and K3 = C, where C is 
nonrandom and depends on the constant K in (j3.9p only. 

Example 3.2 Let Xt g{t,B{t)), where _B is a standard Brownian motion and g{t,x) is a (jointly) 
continuous function on [0,1] x M, having continuous partial derivatives gt(t,x) :— dg(t,x)/dt, gx{t,x) := 
dg{t,x)/dx, gxx{t,x) = d'^g{t,x)/dx^ . By Ito's lemma, 

dXt = g.At,B{t))dB{t) + (^gt{t,B{t)) + ^gxx{t,B{t))^ dt, 

so that X admits the representation p.l|) with at = gx{t, B{t)), bt ~ gt{t, B{t)) + ■^gxxit, B(t)) and the same 
filtration as in the previous example. Assume that 

\gx{t,x)\>Ki, \gx{s,y)-gx{t,x)\<K{\s^t\^/^ + \y-x\), 

for all (t, x), (s, y) £ [0, 1] x K) and some constants < Ki,K < 00. Then X satisfies the conditions in (|3.3p . 



4. Gaussian processes 
4.1. Assumptions 

Let X = {Xt, t E [0, 1]) be a Gaussian process, with zero mean. Without loss of generality, assume Xq — 0. 
Define ap^{k), the variance of A^'"Ar, and Pp,n{k), the correlation coefficient between A^'"X and A^'^^X, 
i.e. 

Let Bh = {Bnit), t G M) be a fractional Brownian motion (fBm) with parameter < < 1, i.e., a Gaussian 
process with zero mean and covariance such that Ei?// (s)i?/f (t) = ^ (I^P'^ + ^ l'^ ^ ^P'^)- Its p— th 
order increments (A^i?//,j G Z) form a stationary Gaussian process, for any p > 1. In particular, the 
covariance function of AjBh = ^]Bh = BhU + 1) - BhU) and A'^Bh = BhU + 2) - 2BH{j + 1) + BhU) 
can be explicitly calculated: 

E[AoBhA,Bh] = 2-i(|.7 + lp^ + |j-l|2^-2bf^), (4.2) 

1^[AIBhA'jBh] - 2-'{-\j+2\^"+Mj + ir-6\.ir + Mj-ir-\j-2\'"). (4.3) 
From Taylor expansion, 

E[AoBhA,Bh] ^ 2H{2H-l)f"-^, 



Y^[AIBhA]Bh\ ^ 2H{2H -1){2H -2){2H -3)j 



2H-A 



as J — >■ cx), and therefore the first increment, (AjBu), has a summable covariance if and only if < < 3/4, 
while the second increment, (A^Bh), has a summable covariance for any < H < 1. 



Introduce the following conditions: 
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(A.l) There exist continuous functions H{t) £ (0, 1) and c{t) > for t e [0, 1] such that Vj G N* 

\ 2 



hm sup 

n->oOtg(0,l) 



{[nt]+j)/n - X[nt]/n 



-c(t) 



1, 



= 0, with 



hm sup \H{t) - H{t + -)\\ogn ^ 0. 

"^°°tG(0,l) " 

(A.l)' There exist continuous functions H{t) e (0, 1) and c{t) > for t e [0, 1] such that Vj G N* 

\2 



hm sup ^/n 

"->~tG(0,l) 



^[^([nt]+j)/n — X[nt]/n) 



2H{t) 



~ 0, with 



(4.4) 
(4.5) 

(4.6) 



hm sup \H{t) - H{t +-)\^\ogn = Q and hm sup \c{t) - c{t + -)\^ ^ Q . (4.7) 
(A.2)p There exist d > 0, 7 > 1/2 and <9 < 7/2 such that for any I < k < j < n with n eW 

E[AP'"^ AP'"X]| < d ap,nik)ap,n{j) ■ n' ■ \j - k\-\ (4.8) 

A straightforward apphcation of Assumption (A.l) (or (A.l)') imphes that \/c{t)BH(t) is the tangent 
process of X for aU t £ (0, 1) and more precisely: 

Property 4.1 Assumptions (A.l), (A.l)' respectively imply that, for any j £ Z and p — 1,2, 



hm sup 



hm -y/n sup 
te(o,i) 



E[Ag^"i?^(,) Af"i3^(,)] 



-c(i) 
-c(t) 



= 0, 



(4.9) 
(4.10) 



Moreover, for any t £ (0, 1) and p = 1,2 



Assumption (A.l) can be characterized as uniform local self- similarity of (Xt) (the uniformity refers 
to the supremum over t £ (0,1) in (j4.4p ). Note that for X having stationary increments and variogram 
V{t) = EXf, Assumption (A.l) reduces to V{t) ~ ct"^" (c > 0, < < 1). For j = 0,1, relation (gH) 
implies that for any t £ (0, 1), the variance and the (1/n)— lag correlation coefficient of Aj'^"jX satisfy the 
following relations: 

/ 1 \ 2H(t) 

^ c{t)a!{H{t))=c{t)E[{AoBHit))']=c{t)(-) , (4.11) 
Pi(ff(t))= corr(i?^(,)(l),i?^(,)(2)-BH(t)(l)) = 2^^(*)~i-l, (4.12) 



C^2,n(M 
P2,n{[nt] 



c(t)a2^(if(t)) =c(t)E[(Agi3^(,))2] ^c(t)(4-4^(*))(i)'''^*\ 
P2iHit)) = corr(Bff(t)(2) - 25^(4), ^^(t) (3) - 25^(4) (2) + Sff(t)(l)) 

_^2H(t) ^ 22-H"(t)+2 _ 7 

8 - 22^(t)+i 



(4.13) 



(4.14) 



(see (|6.31|) '). Moreover, relations (|4.11|) - (|4.14l) hold uniformly in t £ (0, 1). Condition (|4.5p is a technical 
condition which implies (and is "almost equivalent" to) the continuity of the function t — >■ H{t). Assumption 
(A.l)' is a sharper convergence condition than Assumption (A.l) required for establishing central limit 
theorems. 
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Condition (j4.8p specifies a nonasyniptotic inequality satisfied by tlie correlation of increments A^'"X. 
The particular case of stationary processes allows to better understand this point. Indeed, if {Xt) has 
stationary increments, the covariance of the stationary process (A^'"X, fc e Z) is completely determined by 
the variogram V{t), e.g. 

EK^-X AJ.-.Y] = + V(i^) - 2v(tl)}. (4.15, 

In the "most regular" case, when X = Bh is a fBm and therefore V{t) — t^^ , it is easy to check that 
assumption (A.2)2 holds with 61 = and 7 = 4 - 2i/ > 2 (0 < i? < 1), while (A.2)i with 6* 0, 7 = 2 - 2iJ 
is equivalent to H < 3/4 because of the requirement 7 > 1/2. However, for X = Bh, (A.2)i holds with 
appropriate 6' > in the wider region < H < 7/8, by choosing 6 < 2 — 2H arbitrarily close to 2 — 2H and 
then 7 < 2 — 2H + 9 arbitrary close to 4 — AH. A similar choice of parameters 9 and 7 allows to satisfy 
(A.2)p for more general X with stationary increments and variogram V{t) ~ ct'^" (t^O), under additional 
regularity conditions on V{t) (see below). 

Property 4.2 Let X have stationary increments and variogram V{t) ^ ct^^ {t — > 0), withe > 0, H E (0, 1). 
(i) Assume, in addition, that < H < 7/8 and \V"{t)\ < Ct^'^ {0 < t < 1), for some C > and 
4 - AH > K> 2- 2H, K > 1/2. Then assumption (A.2)i holds. 

(a) Assume, in addition, that \V^'^\t)\ < Ct^'^ {0 < t < 1), for some C>0and8-4i7>K>4 - 2H. 
Then assumption (A. 2)2 holds. 

The following property provides a sufficient condition for (A.2)p in spectral terms, which does not require 
differentiability of the variogram. 

Property 4.3 Let X be a Gaussian process having stationary increments and the spectral representation 
(see for instance Cramer and Leadbetter, 1967) 

Xt = I (e'*« - 1) w^(de), for all t e M, (4.16) 

where W{dx) = W{—dx) is a complex-valued Gaussian white noise with zero mean and variance E |VK(da;)|^ = 
dx and f is a non-negative even function called the spectral density of X such that 

^(lA|en/(C)dC < 00. (4.17) 

Moreover, assume that f is differentiable on (K, 00) and 

fiO-cC^"-^ (C^oo), \f'm<Cr^"-' {^>K) (4.18) 

for some constants c,C,K > 0. Then X satisfies assumption (A.2)i for < H < 3/4 and assumption 
(A.2)2 for 0<H <l. 



4.2. Limit theorems 

Before establishing limit theorems for the statistics RP'" for Gaussian processes, recall that A is given in 
(11.91) and with Pp{H) in (16.311) one has 

/ Ap{H{t))dt = / X{pp{H{t)))dt - / E 
Jo Jo Jo 

Straightforward computations show that Assumptions (A.l) and (A.2)p imply Assumption (A) with 
A(*)((5) = ,5-^(*),y(*) = y^dJ)BH(t)- Therefore Proposition [l^l] ensures the convergence (in L^) of the 
statistics RP'"" to Ap{H{t))dt. Bardet and Surgailis (2009) proved a.s. convergence in Theorem 14.11 
below, using a general moment bound for functions of multivariate Gaussian processes (see Lemma [1] in 
Section [5]). A sketch of this proof can be found in Section [51 



I^oBhh) - 


^ ^lBH(t) 1 




^ \^lBH(t)\ 



dt. 
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Theorem 4.1 Let X be a Gaussian process satisfying Assumptions (A.l) and (A.2)p. Then, 

^ [\p{H{t))dt {p = l,2). (4.19) 



Corollary 4.1 Assume that X is a Gaussian process having stationary increments, whose variogram satisfies 
the conditions of Properties\J7^ or\4.S\ Then 



^P,n ^ (p=l,2). (4.20) 

n— >oo 

The following Theorem 14. 2 1 is also established in Bardet and Surgailis (2009). Its proof (see a sketch of this 
proof in Section [6]) uses a general central limit theorem for Gaussian subordinated nonstationary triangular 
arrays (see Theorem[T]in Section[6]). Note that the Hermite rank of 'iIj{xi,X2) = \xi + a;2|/(|a;i| + |a;2|) is 2 
and this explains the difference between the cases p — I and p — 2 in Theorem 14.21 in the first case, the 
inequalities in (|6.12p for (Yn{k)) as defined in (|6.9p - (l6.10p hold only if suptg[Q H{t) < 3/4, while in the 
latter case these inequalities hold for < sup^^jQ ^ Hit) < 1. A similar fact is true also for the estimators 
based on generalized quadratic variations, see Istas and Lang (1997), Coeurjolly (2001). 

Theorem 4.2 Let X be a Gaussian process satisfying assumptions (A.l)' and (A.2)p, with 6 — 0. More- 
over, assume additionally sup H(t) < 3/4 if p = 1. Then, for p — 1,2, 

tG[04] 

(4.21) 



lUp.n_ f\^(^Hit))dt) A Af(0, f\p{H{T))dT), 
V Jo n^oo ^ Jo ' 



with Ap(H) and Y,p{H) given in il.^) and hi .12]) . respectively. 

The following proposition shows that the previous theorems are satisfied when smooth multiplicative and 
additive trends are considered. 

Proposition 4.1 Let Z = [Zt = a{t)Xt + f3{t), t E [0,1]), where X = {Xt,t E [0,1]) is a zero mean 
Gaussian process and a, (3 are deterministic continuous functions on [0, 1] with infjg[o,i] ctit) > 0- 

(i) Let X satisfy assumptions of Theorem \4-.2\ and a E C[0, 1], /3 E C^[0, 1]. Then the statement of Theorem 
\4.2\ holds with X replaced by Z. 

(ii) Let X satisfy assumptions of Theorem \4.1\ and a E C^[0,1], /S E C^[0, 1]. Then the statement of Theorem 
T] holds with X replaced by Z. 



Remark 4.1 A version of the central limit theorem in (|4.2ip is established in Bardet and Surgailis (2009) 
with Ap{H{t))dt replaced by Ei?^'" under weaker assumption than (A.l)' or even (A.l): only properties 
(|4TT|) - (|4J2|) (for p=l) and (|4l^ - (j4lll) (for p = 2), in addition to (A.2)p with 6* = , are required. 

The particular case of Gaussian processes having stationary increments can also be studied: 

Corollary 4.2 Assume that X is a Gaussian process having .stationary increments and there exist c > 0, 
C > and < H < 1 such that at least one of the two following conditions (a), (b) hold: 

(a) variogram V{t) ^ ct^" {l + o{t^l'^)) for t and \v'^^p'> {t)\ < Ct^^-^P for all t E (0,1]; 

(b) .spectral density f .satisfies (4H\ ), (4l8\ ) and f{Q = c(-^"-^{l + 0(^-^/2)) (^ ^ 00). 

Then: 

MRP^--ApiH)) ^ M(o,MH)) = (4-22) 

ri^oo \ / IP — 2, < H < 1. 

Moreover, with the expression of s\[H) given in Section\^ 

Mk-\R^-^^)-h) a N{Q,sl{H)). (4.23) 
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Therefore, -ff„ — ^^'^{R^'^) is an estimator of the parameter H following a central limit theorem with 
a convergence rate ^Jn under semi-parametric assumptions. Similar results were obtained by Guyon and 
Leon (1989) and Istas and Lang (1997) for generalized quadratic variations under less general assumptions. 

Remark 4.2 In the context of Corollary 14.21 and H E (3/4, 1), we expect that R^'" follows a nongaussian 
limit distribution with convergence rate n^^^^ . 

Remark 4.3 Figure 1 exhibits that H i— )■ A2{H) is nearly linear and is well- approximated by .1468if-f .5174. 
Consequently, A2(ff(t))dt w .1468//+ .5174, where H = H{t)dt is the mean value of the function H{-). 



Another interesting particular case of Theorem 14 . 2 1 leads to a punctual estimator of the function H(t) from 
a localization of the statistic R^'". For tg S (0, 1) and a E (0, 1), define 



[nto+r, 



1 



This estimator is studied in Bardet and Surgailis (2010) and compared to the estimator based on generalized 
quadratic variations discussed in Benassi et al. (1998) and CoeurjoUy (2005). 



4.3. Examples 

Below we provide some concrete examples of Gaussian processes which admit a fBm as the tangent process. 
For some examples the hypotheses of Theorems 14 . Hl¥^ and the subsequent corollaries are satisfied. For other 
examples, the verification of our hypotheses (in particular, of the crucial covariance bound (A.2)p) remains 
an open problem and will be discussed elsewhere. 

Example 4.1 Fractional Brownian motion (fBm). As noted above, a fBm X = Eh satisfies (A.l)' as well 
as (A.2)i (for < i/ < 3/4 if 61 = and < iJ < 7/8 if < 6* < 2 - 27? with 9 arbitrary close to 2 - 2iJ 
and therefore j < 2 - 2H + 9 arbitrary close to 4 - AH may satisfy 7 > 1/2) and (A. 2)2 (for < H < 1), 
with H{t) = H, c{t) = c. Therefore, for fBm both Theorems 14.11 (the almost sure convergence, satisfied for 
< H < 7 /8 when p — 1 and for < H < 1 when p = 2) and 14.21 (the central limit theorem, satisfied 
for < < 3/4 when p = 1 and for < < 1 when p — 2)) apply. Obviously, a fBm also satisfies the 
conditions of Corollary 221 Thus, the rate of convergence of the estimator A^^iJ^.n) ^. jj^ jj 

is \/n. But 

in such a case the self-similarity property of fBm allows to use in this case asymptotically efficient Whittle 
or maximum likelihood estimators (see Fox and Taqqu, 1987, or Dahlhaus, 1989). However, for a fBm with 
a continuously differentiable multiplicative and additive trends, which leads to a semi-parametric context, 
the convergence rate of i7„ is still ^/n while parametric estimators cannot be applied. 

Example 4.2 Multiscale fractional Brownian motion (see Bardet and Bertrand, 2007) defined as follows: for 
£ G N*, a (M£)-multiscale fractional Brownian motion X — {Xt, t eM.) {{Mi)- fBm for short) is a Gaussian 
process having stationary increments and a spectral density / such that 

2 

/(^) = ^5j+T IK ^ 1^1 < ^^+1) fo"^ all ^ e M (4.24) 

with luq := < uji < • • • < LOg < w^+i :— 00, ai > and Hi E M for i E {0, • • • , ^} with Hq < 1 and 
Hi > 0. Therefore condition (|4.18p of Propertv 14.31 is satisfied, with K = uj( and H = H(. Moreover, the 
condition /(^) = c^^^^^^(l + o(f^^/^)) (^ -> 00) required in CoroUarv 14.21 is also checked with H = He. 
Consequently, the same conclusions as in the previous example apply for this process as well, in the respective 
regions determined by the parameter Hg at high frequencies x > ujg alone. The same result is also obtained 
for a more general process defined by /(^) = c^~^^~^ for |^| > ui and condition (I4.17|) is only required 
elsewhere. Once again, such conclusions hold also in case of continuously differentiable multiplicative and 
additive trends. 
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Example 4.3 Multifractional Brownian motion (mBm) (see Ayache et al., 2005). A mBm X — {Xt, t € 
[0, 1]) is a Gaussian process defined by 

Xt - = g{H{t)) J^j^^^Widx), (4.25) 

where W^(da;) is the same as in (j4.16p . H{t) is a (continuous) function on [0, 1] taking values in (0, 1) and 
finally, g{H{t)) is a normalization such that EXf = 1. It is well-known that a mBm is locally asymptotically 
self-similar at each point t € (0, 1) having a fBm i?ff(t) as its tangent process at t (see Benassi et al, 1997). 
This example is studied in more detail in Bardet and Surgailis (2010). 

Example 4.4 Time-varying fractionally integrated processes. Philippe et al. (2006, 2008) introduced two 
classes of mutually inverse time- varying fractionally integrated filters with discrete time and studied long- 
memory properties of the corresponding filtered white noise processes. Surgailis (2008) extended these filters 
to continuous time and defined "multifractional" Gaussian processes {Xt, t>Q) and {Yt,t > 0) as follows 



Xt 



T{H{1)-.^) ^' .s)^^(-'-i-^e^-(-^)dr}dB(.), (4.26) 
- .)f ^^-^e-^-(-*) - (-.)f ^)-^e-^+(-°)}di?(.), (4.27) 



where s" := s"l(s > 0), -B is a Brownian motion, 

, , , /■* H{u) - H{t) , . , ^ /"* H{s) - H{v) , 
A^{s,t):= / !^dM, A+{s,t):^ / ^-^dv {s < t) 



and where H{t),t G R is a general function taking values in (0, oo) and satisfying some weak additional 
conditions. Surgailis (2008) studied small and large scale limits of {Xt) and {Yt) and showed that these 
processes resemble a fBm with Hurst parameter H = H{t) at each point i e K (i.e., admit a fBm as a 
tangent process) similarly to the mBm in the previous example. The last paper also argues that these 
processes present a more natural generalization of fBm than the mBm and have nicer dependence properties 
of increments. We expect that the assumptions (A.l), (A.l)', (A.2)p can be verified for (|4.26|) . (|4.27|) : 
however, this question requires further work. 



5. Processes with independent increments 

In this section, we assume that X = {Xt, t > 0) is a (homogeneous) Levy process, with a.s. right continuous 
trajectories, Xq — 0. It is well-known that if the generating triplet of X satisfies certain conditions (in 
particularly, if the Levy measure v behaves regularly at the origin with index a G (0,2)), then X has a 
tangent process Y which is a—stable Levy process. A natural question is to estimate the parameter a with 
the help of the introduced statistics W'^ . Unfortunately, the limit of these statistics as defined in (1.5) 
through the tangent process depends also on the skewness parameter /3 e [—1,1] of the a— stable tangent 
process Y and so this limit cannot be used for determining of a if /3 is unknown. 

In order to avoid this difficulty, we shall slightly modified our ratio statistic, as follow. Observe first 
that the second differences A^'"A' of Levy process have a symmetric distribution (in contrast to the first 
differences A^'"X which are not necessary symmetric). For notational simplicity we shall assume in this 
section that n is even. The modified statistic 

is written in terms of "disjoint" (independent) second order increments {H^^ X, AjjI^^j"'^) having a symmetric 
joint distribution. Instead of extending general result of Proposition 2.2 to i?^'", we shall directly obtain its 
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convergence under suitable assumptions on X . Note first 

Ei?^- = E^(xj5) -X(Ji), (5.1) 

where X^^\i = 1, • • • , 4 are independent copies of X. Note tliat 1/2 < Ei?^'" < 1 since 
Ei;(X<:^^ ^X\]KX<:^;> - xj^M > P(xj^) _ X « >0,xj^) - X^^^ > 0) 

1/n 1/n' 1/n 1/n/ — ^ 1/n 1/n — ' 1/n 1/n — / 

+ p(x|^) <o,x(^) -xj^) <0) 

y 1/n 1/n ' 1/n 1/n / 

> P'ix[% - Xjji > 0) + P2(X(% - X[l < 0) > 1/2. 



1/ 

Proposition 5.1 Let there exists a limit 



The 



lim E^2,n ^ (5 2) 

n— )-oo 



n— >oo 



Proof. Write i?^'" = E7?^'" + (n/2 — 1)^^(3ti, where (3„ is a sum of centered 1-dependent r.v.'s which 
are bounded by 1 in absolute value. Therefore E((n/2 — l)^^Qn)^ — 0{n^^) and the a.s. convergence 
{n/2 — l)^^Q,i — > follows by the Chebyshev inequality. □ 

Next we discuss conditions on X for the convergence in (|5.2p . Recall that the distribution of is infinitely 
divisible and its characteristic function is given by 



Ee'^^' = exp 



|t(i76l-ia2e'2+y" (e™^-l-iu6ll(|u| < l))i/(du))}, 61 e M, (5.4) 



where 7 G M, a > and v is & measure on R such that Jj^ min(u^, l)z^(du) < 00. The triplet {a,^,v) 
is called the generating triplet of X (Sato (1999)). Let X'^'^\i = 1,2 be independent copies of X. Note 
is a Levy process having the characteristic function 



Ee'^^' =exp|i(-a26)2 + 2^ Re(l - e™'')dX(u)) }, e M, 



(5.5) 



where 

K{u) :— j/((— cx), — u] U [u, 00)) 

is monotone nonincreasing on (0, 00). Introduce the following condition: there exist < a < 2 and c > 
such that 

Kiu)^^, uiO. (5.6) 

It is clear that if such number a exists then a :— inf {r>0; |a;|'':^((ix) < 00} is the so-called fractional 

order or the Blumcnthal-Getoor index of the Levy process X. 

Let Za be a standard a— stable r.v. with characteristic function Ee'^^° = c^'^' and Za\i = 1,2,3 be 
independent copies of Za- 

Proposition 5.2 Assume either a > or else, a — and condition 115. 6\) with < a < 2 and c > 0. Then 

^-i/q(^U) _ Xj2^) cZa with c depending on c, and \5.2^ . \5. g)] hold with 
t-yo 

A^A(a) :=E^(ZW,Z(2)). 
Moreover, with a^{a) := 2vaT{4>{Zi^\ Z^J^)) + 4cov{^p{zi^\ Z^^^) ,^{zi^\ Zi^'^)) , 

V^(ii2.» - Eii^'") A NiO,a^{a)). (5.7) 
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Proof. Relation t ^^"Wt — t ^/"(Xj — ) — >■ cZa is an easy consequence of the assumptions of the 

proposition and the general criterion of weak convergence of infinitely divisible distributions in Sato (1991, 
Theorem 8.7). It implies (j5.2[) by the fact that is a.e. continuous on M^. Since i?^'" is a sum of 1-dependent 
stationary and bounded r.v.'s, the central limit theorem in (j5.7l) follows from convergence of the variance: 

nv£iT{R^-") ^ a^{a); (5.8) 

see e.g. Berk (1973). Rewrite R^^"- = (n/2 - 1)-^ J^klTa^^^^ Vn{k), ?7„(fc) := ip{Al'^X, AI'^I^X). We have 

nvar = var 77,, + \' cov r?„ , 7?„ 1 , 

n/ 2 — 1 (."■/2 — 1) 



where var(7)„(0)) 
as in the proof of 



var(^(zi'\zi'))), cov(77„(0),77„(l)) cov(V'(zi'\ zi'^) , ^i^^)) similarly 

n—^oo 

above. This proves f|5.8p and the proposition. □ 




0.4 0.6 



Figure 4. The graphs of a 1— > A(q^) = E ' ,°. %J- (left) and a 1— > cr(Q!) (right) for a process with independent increments. 



The graph of A(a) is given in Figure H) Note that A(2) — Ai(l/2) ~ 0.72: this is the case of Brownian 
motion. 

In order to evaluate the decay rate of the bias Ei?^'" — A{a) we need a uniform convergence rate in Lemma 
15.11 below, for 

\\Fn ~ Galloo sup \Fnix) - G„(a;)|, F„(x) := P(ni/"I^i/„ < x), G„(a:) := P(Z„ < x). 



where := cZq, is the limiting a— stable r.v. in Proposition 15.21 and {Wt,t > 0) is the symmetric Levy 
process with characteristic function as in (15. 5p . The proof of Lemma |5. II is given in Annexe. 



Lemma 5.1 (i) Let a = and K satisfy 15. 6]) . Denote Ki{u) :— K{u) — cu ji^TiKu) := \(lKi{v)\, 
the variation of Ki on [u^oo). Moreover, assume that there exist some constants /3,(5 > such that 

|ifi|(ii) = 0(ii-("-'^)+) (ii->0), \Ki\{u) ^ 0{u-^) (ii->c5o), (5.9) 

where x+ :=max(0, x). Then 

{0(n-^/"), ifP<a, 
O(n-ilogn), ifP = a, (5.10) 
O(n-i), ifl3>a. 

(ii) Let a > and K satisfy 

K{u)=0{u-°') (u^O), K{u) = 0{u-^) (w^oo), (5.11) 
for some < a < 2, 5 > 0. T/ien 

1^ (^(n logn), ija — 0. 
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Proposition 5.3 Assume either a > or else a = and condition 115.6]) . Then for any a € (0,2] 

/•oo 

lEii^."- A(a)| < 2C||F„-Ga|loo, ^ := {1 + z)-^dz < oo. (5.13) 

Proof. Let ipi^, y) :— \x — y\/ {x-\-y), x,y > 0^ and let F„, Ga be the same as in Lemma [5?T] Similarly as in 
Vaiciulis (2009, proof of Th. 1), write 



OO pOO 



Ei?2."_A(a) = 2/ / i;{x,y){dFn{x)dFn{y)-dGa.{x)dGo.(y)) = 2{Wi + W2), 
Jo Jq 

where Wi := J^"" ^^(x, y)di^^„(a;)(dF„(2/) - dG„(y)), W2 := V^(:c, y)dG„(2/)(dF„(x) - dG„(x)). 

Integrating by parts yields 

iH^il = 2 r \x\dF^{x) r \F^{y) - G^{y)\^^^ 
Jo Jo + V) 

dy_ 

[x + yy 



< 2\\Fr,-G^\\oc I \x\dFn{x) I ' = GHf„-G^Hc 



since dF^ix) = 1/2. A similar estimate holds for W2. This proves (|5.13p . □ 

Propositions 15. 2[ 15.31 and Lemma [5.11 together with the Delta- method, yield the following corollary. 

Corollary 5.1 Let a and K satisfy either the assumptions of Lemma 5.1(i) with /? > a/2, or the assumptions 
of Lemma 5.1 (ii). Then 

V^(^2'" - A(a)) A Af{0,a\a)). 



Moreover, if we define (in := A ^(K^'"), then 

y/n{an-a) 7V(0,S^(a)), 



n— J-oo 



where s^{a) := [f^(a)] a^{a), 0<a<2. 

There exist very few papers concerning estimation of a in such a semiparametric frame. Nonparametric 
estimation of parameters of Levy processes based on the empirical characteristic function has recently been 
considered in Neumann and Reifi (2009) and Gugushvili (2008), but the convergence rates there are (logn)'' 
with K > 0. Ait Sahalia and Jacod (2009) have proposed an estimator of the degree of activity of jumps (which 
is identical to the fractional order in the case of a Levy process) in a general semimartingale framework using 
small increments of high frequency data. However from the generality of their model, the convergence rate of 
the estimator is not rate efficient (in fact smaller than n^/^). A recent paper of Belomestny (2010) provides 
an efhcient data-driven procedure to estimate a using a spectral approach but in a different semiparametric 
frame from ours. Thus, CoroUarv 15.11 appears as a new and interesting result since the estimator a„ follows 
a -v/n-central limit theorem. 



6. Annexe : proofs 

Proof of Lemma \3.1\ 

Let 6^ := maxj=i,2 If > 1/2, then holds since the l.h.s. of ([32]) does not exceed 1. Let (5^ < 1/2 
in the sequel. Write U V(^i + 6, ^2 + 6) - V'(^i, ^2) ^ Us + Ul, 

Us Ul{As) - (V'(Zi+ei,^2+6)-V'(^i,^2))l(^5), 

Ul U\{A<i) = (V(Zi+a,^2+6)-V'(^i,^2))l(A^), 
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where 1(^5) is the indicator of the event 

As {\Zi\ > 52/3, IZ2I > 5^/', 161 < ^'/V2, 161 < 5^/^/2}, 
and 1(^5) = 1 — 1(^5) is the indicator of the complementary event A"^. Clearly, 

E|C/|| < 2(P(|Zi| < (5^/3) + P{\Z2\ < <52/3) + P(|6| > 6^/y2) + P(|6| > 6^/y2) 

< -^Sy^ + 2^ <4j2/3. 
~ V2^ (54/3 - 



It remains to estimate EjL'al. By the mean value theorem, 



lUsl < la I sup I V-. J + 16 1 sup 1^, J 1(^5), 

D D 



where 



supji/'ajj := snp{\d'ip{xi,X2)/dxi\ : \xi - Zi\ < ^ 1,2}, 

D 

SUplV'ajal 
D 



snp{\dtlj{xi,X2)/dx2\ : \xi - Zi\ < |6|,i = 1,2}. 



Therefore 



Next, 



Therefore 



< s{E'/'[{snp\^,,\)H{As)] +EI/2 [(sup 1^,, 1)21(^5)] 

D D 



\sgn{xi + X2){\xi\ + \X2\) - i\xi + X2\)sg-n{x^)\ 2 

\iPxi(Xl,X2)\ = -. : : < 



(l^l| + l^2|) = 



a;i + \x2\ 



snpl^P^fliAs) < 4sup{(|xi| + |x2|)-2 : \x, ^ Z,\ < d^/y2,\Z,\ > S^/^i ^ 1,2} 



D 



< 16(|Zi| + |Z2|)-2l(|Z,| > ,52/3,1 = 1,2), 



implying 



where 



E[(sup|V'.J)'l(^5)] < 16E 

D 



{\Zl\ + \Z2\f 



■\Z,\>5^/\i^ 1,2 



< C{S), 



CiS) < ^ / ^-^e-(-?+-i)/2da;idx2 
27r J{xj+xl>S'i/^} ^1 + ^2 



< 16 / r-^e-'^'/^dr < 16 (1 + (4/3)| log 

Hence E\Us\ < 8 (5(1 + (4/3) | log(5|)i/2, E\U^\ < 4(5^/3. It remains to use x{l + (4/3) | loga;|)i/2 < 2x2/3 for 
all < X < 1. □. 



Proof of Property 4 ■ 1 

We use the following identity: for any reals xi, - ■ ■ ,Xj, 



fc=l fc=2 k=l 



^ X 2 
fe=2 
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In particular, 



E 



[nt] 3 + [nt] 



- 2 1 



l+[nt.] 



where :— t + (so that [nt,] = [nt] + 1). Then, using (A.l) and the notation u„ for a sequence 

tending to as n — )■ oo uniformly in t and all \ j\ < J, where J is a fixed number, we obtain 



E 



1 r / \2Hlt) / \2Hlt,) 

/ .\ 2H(t) / .X 2H(t.) ^ 



/ .N 2ff(t.) / .\2H{t) 

since c(t*) — c(t) = w„ and f ^ j ~ ( n ) (-'^ '^") foUo'^s from (|4.5p . This proves (|4.9p for p = 1. 

Relation (|4.9p for p = 2 follows analogously. Relation (|4.10p also follows by the same argument and the 



fact that c(i*) — c{t) ~ Unj \/n and 
Property 14. II is proved. 



2H{t, 



/ .\2ff(t) 

I ^ j (1 + Un/ \/n) hold in view of Assumption (A.l)'. 



□ 



Proof of Property \4.S\ 

With condition V{t) ^ ct^^ {t — >■ 0) in mind, inequality (|4.8p reduces to 

\V (^ii) + V (^) - 2V (^) I < Cn-^^+^/fc-^ (p = 1, 2 < fc < n), (6.1) 

\V {!^) -AV {^) +6V {^) -W {!^) +V {!^)\ < Cn~^"+'^k-^ (p = 2, 4 < fc < n). (6.2) 

The left hand side of (16.11) can be written and estimated as 



l/n f\l 



JQ 



V"{t - s^{kln))AtAs 



/•I/" /■!/" 

< C / / |t-s + (fc/n)|-''dMs 



1 /.I 



< Cn« 



|t-s + /s|-'"dMs 



-2//i,-7 



where 7 = k > 1/2 and 9 = k + 2H — 2 G [0, 7/2) since k < 4 — 4_ff. This proves part (i). Part (ii) follows 
similarly, by writing the left hand side of (|6.2p as 



l/n /.1/ri 



Jo 



V'-'^\t -S + U-V + {k/n))dtdsdudv 



/.l/n /.l/n 

< C ■ ■■ \t- s + u~ V + {k/n)\-'^dtdsdudv 

Jq Jo 



< Cn^-^k-" = Cn'^-^" k-'-' , 
where 7 = k > 1 and 6 = k + 2H - 4 e [0, k/2) since k < 8 - 4i/. Property 112] is proved. 



□ 



Proof of Property 
From (|4.16l) we have 



,i(2;/n) 



2p . 



\x{j/n 



^f{x)dx = 2^+p / (1 - cos(x/n))Pcos(a;j7n)/(a;)dx 
Jo 



2i+f(/i+/2 
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\h\ 



K 



(1 — cos(a;/n))^ cos{xj/n)f{x)dx 



C 





Cn-^^ [ 




Jo 



K 



(6.3) 



and h = n (1 — cos(x))p cos{xj)f{nx)dx 

JK/n 



(n/j) / (1 - cos(a;))^'/(7ia;)dsin(a;j) 

'K/n 



= -{n/j) / sm{x3){{l~cos{x))Pf{nx))'^dx + 0{n-^P) (6.4) 

JK/n 

since \{n/j)f{K){l - coii{K / n)Y K / n)\ < C{n/j){K/nYP\jK/n\ < Cn-^P for any ivT > fixed. 
Let p — I- The last integral can be rewritten as sin(a;j) ((1 — cos{x)) f (nx))'^ dx — Ii + I2, where 

/>oo poo 3 

< / I sin(xj) sin(a;)/(na;)|da; < C / | sin(xj) sin(a;)|(na;)^^^^^da; < ^^^Iiq, 

JK/n JK/n i 



g=l 



where we used the fact that f{x) < Cx ^ {x > K); see condition (|4.18p . and where 













Jo 


/12I 




/"' 

poo 


/13I 




I 1 



\sm{xj)sm{x)\{nxy'^"-^dx < Cj I x^{nx)-'-'^" dx < Cf"-'n 



■2H-1„-1-2H 



|sin(xj)sin(x)|(nx)-^-""Ma; < C / x{nx)-'-'"dx < Cj^"-'n 



1/j 



|sin(xj)sin(x)|(nx)"2^"Ma; < C / {nxy-~'^" dx < Cn 



-1-2H 



Similarly, using (|4.18p . 

poo POO 3 

I/2I < n |sin(a;j)(l -cos(x))/'(na;)|da; < C | sin(a;j)(l - cos(x)|(nx)"^^"2dx < V/29, 

JK/n JK/n i 



where 



I21 = n \sm{xj){l-cos{x))\{nx)-^"-''dx < Cnj x\nx)-^-^^dx < Cj^"-^n 



q=l 



:2H-l„-l-2H 



I22 = n \sm{xj){l-cos{x))\{nx)-''"-Mx < Cn I x\nx)-^~^"dx < Cj^^-'^n 



2H-1-1-2H 



1/j 

00 



J23 = y I sin(a;j)(l - cos(a;))|(na;)"^^"Ma; < Cn / (na;)"^^^^da; < Cn 
We finally obtain, for 1 < < n. 



-1-2H 



2H-2 



implying (for < iJ < 3/4) (A.2)i with 6 = and k = 2 - 2H > 1/2. For p = 2, the estimation of the 
integral in ()6.4p is completely similar, resulting in the bound 



for any < < 1, or (A.2)2 with 6* = and k = 1 > 1/2. 



□ 
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Sketch of the proof of Theorem \4.1 



The proof of Theorem 14.11 is based on the moment inequality in Lemma [H below, which extends a similar 
inequality in (Taqqu, 1977, Lemma 4.5) to vector-valued nonstationary Gaussian processes. The proof of 
Lemma [T] uses the diagram formula and is given in Bardet and Surgailis (2009). To formulate this lemma, 
we need the following definitions. Let X be a standard Gaussian vector in R'^ {ly > 1) and let L^(X) 
denote the Hilbert space of measurable functions / : M'^ M satisfying :— E(/(X))^ < oo. Let 

L§(X) = {/ G L^(X) : E/(X) = 0}. Let (Xi, • • • ,Xn) be a collection of standardized Gaussian vectors 
Xt = {xj:^\ ■ ■ ■ jXj:"'') G M.'^ having a joint Gaussian distribution in R''^. Let e € [0, 1] be a fixed number. 
Following Taqqu (1977), we cah (Xi, • • • , Xn) e-standard\i \Ex[''^x'f^\ < e for any t ^ s,l < t, s < N and 
any I < u,v < i^. Finally, denotes the sum over all distinct integers I < ti, ■ ■ ■ ,tp < N, ti ^tj {i =/= j). 

Lemma 1 Let (Xi, • • • ,Xjv) be e—standard Gaussian vector, Xt ~ {x[^\ ■ ■ ■ jXj"^) e R"^ [v > 1), and let 
Gj,t,N € L^(X), 1 < J < p (p > 2), 1 < t < A^. For given integers m, N > 1, define 

.^f^A E ,"1^^, |Ex(")xi")r. (6.5) 

l<t<N ^ — ^ l<u,v<i' 
l<s<N,s=jLt ^ 

Assume that for some integer Q < a < p, the functions Gi^t,N i ' ' ' i Ga^t.N have a Hermite rank at least equal 
to m for any N > 1,1 < t < N , and that e < ^ . Then 



^'E|Gi,t„Ar(XtJ...Gp,t^,^(XtJ| < C[e,p,m,o,,v)KN^- 
where the constant C{s,p,m, a, v) depends on e,p,m,a,h' only, and K := Y[''j^iT^SLKi<t<N \\Gj^t,N\ 



Sketch of the proof of Theorem \4-l\ The convergence lim„^oo Ei?^^" = Ap{H {t))dt is easy (see the proof 
of Proposition 2.2). Hence (|4.19p follows from 

:= - Ei?P^" ^ 0. (6.6) 

Relation (|6.6p follows from the Chebyshev Inequality and the following bound: there exist C,k > 1 such 
that for any n > 1 

e(rP'A < Cn-". (6.7) 



By definition, = ^ EVo^' ^n{k), where 77„(fc) := 77„(fc) - Er,„(fc) and ?7„(fc) := V(A^'"^, AP:;^iX), 
tp{x,y) — \x + y\/{\x\ + \y\) are nonlinear functions of Gaussian vectors (A^'"X, A^^^j^X) G R^ having the 
Hermite rank 2; however, these vectors are not e— standard and therefore Lemma 1 cannot be directly applied 
to estimate the l.h.s. of ()6.7p (with p = 1, ■ ■ ■ ,A, v = 2). To this end, we first need to "decimate" the sum 
RP'^, as follows. (A similar "trick" was used in Csorgo and Mielniczuk (1996).) Let £ = [n^^'^] be the 
sequence of integers increasing to oo (at a rate o(n^/^) by condition 6 < 7/2) and write 

£-1 . [{n-2-])/e] 

j=0 k=0 

Then 

4 



e(^-) <.^maxE(Er(.)) 



0<j<£ 

Write rjnik) as a (bounded) function in standardized Gaussian variables: 

rin{k) = fk,n {Yn{k)) , 
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where Y^{k) = (ri''(fc),yi')(fc)) e 

Yjl'Hk) (6.9) 
fs-w X ^V'X p„Jk) Al'"X 1 ^ ^ 

and A-,„(x«,x(^)) ^ (^^'^ '''t^)^ [PvAk)x^'^ + ' PlAk)^^'^)) ■ (S-H) 

where (Tp „(fc), Pp,n{k) are defined in (|4.ip . Then, for each fc, l^n(fc) {Yn^\k),Yn'^\k)) has a standard 
Gaussian distribution in and ?7„(fc) = fk,n {Yn{k)) — Efk,n {Yn{k)). Moreover, the vector (Y'„(H+j), fc = 
0, 1, • • • , [(ri — 2 — j)/£]) e IR2([("-2-i)/f]+i) jg £_standard provided £ is large enough. Now Lemma [1] can be 
used and it implies the bound (|6.7p using Assumptions (A.l) and (A.2)p. The details of this proof can be 
found in Bardet et Surgailis (2009). □ 

Sketch of the proof of Theorem 



The proof of Theorem 14.11 uses the following central limit theorem for Gaussian subordinated multidimen- 
sional triangular arrays. Theorem [T] is proved in Bardet and Surgailis (2009). It extends the earlier results 
in Breuer and Major (1983) and Arcones (1994). Below, similarly as in Lemma 1^ X £ W designates a 
standard Gaussian vector. 

Theorem 1 Let (^n(fc))i<fe<„ be a triangular array of standardized Gaussian vectors with values in 
Yn{k) = (yi^'(fc),---,ri''^(fc)), Eri^'(fc) = O, YYi^\k)Y,[''\k) = 5^^. For a given integer m > 1, 
introduce the following assumption: there exists a function p ; N — > M such that for any 1 < p^q < v , 

V(j,fc)e{l,...,n}2, |EriP)(j)r„(^)(fc)| < \p{j-k)\ with ^|p(j)|™<oo. (6.12) 

Moreover, assume that for any t G [0, 1] and any J G N*, 

(F„([nr]+j))_^<^.<^ A {Wrm-j<,<j, (6.13) 

— n— s-oo — 

where (WV(j))jg2 is a stationary Gaussian process taking values inW and depending on parameter t G (0, 1). 
Let fk,n G L§(JC) {n > 1,1 < k < n) be a triangular array of functions all having Hermite rank at least m. 
Assume that there exists a ]Lq{X)— valued continuous function <Pt,t G [0, 1], such that 

sup ||/V„],„-0.|r = sup E{f\rn].n{X) - MX)f ^ 0. (6.14) 
rG[0,l] tG[0,1] n^oo 



Then, with ^ j dr(^E[0r {Wr{Q))4>r {Wr{j)) 



< oo. 



n-'/'V/fc,„(r„(fc)) A AA(0,a2). (6.15) 

— n— ^oo 

fc=l 

Sketch of the proof of Theorem \4.2\ It suffices to show that 

y^Ei?"'f- / kj,{H{t))dt (6.16) 

and 

v^(i?P^" - Ei?P^") A A/'(0, / Ep(i/(T))dr). (6.17) 
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The proof of (|6.16l) uses Assumption (A.l)' or (I4.10p and the easy fact that for Gaussian vectors {Zn \ Z^^) S 
n e N with zero mean EziP = 0, i 1, 2, n e N and E(Z^^^)2 = E(Z^^')2 = 1, |EZ^^^Z^^^| < 1 

|E^(Z«zf )-E^(zW,Z(2))| < |EZ«Z(^'-ZWZ«|. (6.18) 

The proof of (|6.17p is deduced from Theorem [1] with the sequence of standardized Gaussian vectors Yn{k) = 
{Yj^\k),Yji'^\k)) {v = 2) given in (|6J]) - ((O0| and the centered functions 

/fe,„(a;«,x(2)) :=/fc,„(a;«,a;(2))_E/fe,„(r4fc)), M^^^K^^^^) := 0r(a;W, a;^^)) _ E0,(X) 

with fk.n : — > M given in (|6.1ip and the (hmit) function 

Thanks to symmetry properties of these functions, it is clear that the Hermite rank of fk^n (for any k and n) 
and 4>T (for any r e [0, 1]) is m = 2. Using Assumptions (A.l)' and (A.2)p (with 6* = 0,7 > 1/2), one can 
show that the conditions of Theorem [1] are satisfied for the above fn.k, 4>t and the hmit process [W r{i))j£Z 
in (j6.13p is written in terms of increments of fBm (_B^(^)(j'))jgx: 

Wr{j) := (APi3H(r)(j), (-Pp(i/(r))A?BH(.)(j-) + A?Bh(.)(j- + 1)) 1^1 - pI{H{t))) 

having standardized uncorrelated components. The details of this proof can be found in Bardet and Surgailis 
(2009). □ 



Proof of Proposition 4 ■ 1 



(i) Denote Zt := a{t)Xt, R"'P := i?"^P(Z), ^"'P := E^^p{Z). Clearly, part (i) (the CLT in (li^TIl for trended 
process Z = aX + (3) follows from the following relations: 

pi 

^{RP^'' ^ERP'"") A Af{0, I]p(i/(T))dT), (6.19) 

n-s-oo Jq 



E^P'" - / Ap{H{t))dt 0, (6.20) 

Jo n^oo 

^^|Ei^"■''-E7?"^^'| — > 0, (6.21) 

n— )-oo 

nvar(i?"'''-^"'P) 0. (6.22) 

The central limit theorem in (|6.19p follows from Theorem 221 and Remark |4 . 1 1 since Assumption (A.2)p and 
the convergences (|4.1ip - ()4.14p (with c{t) replaced by a'^{t)c{t)) can be easily verified for the process Z = aX. 

Let us turn to the proof of (I6.19l) - (|6.22l) . For concreteness, let p = 2 in the rest of the proof. Since 
A^'"Z = a(fc/n)A^;"X + 2a'{k/n)n-^Al'^^X + 0{l/n'^)X{{k + 2)/n), it follows easily from Assumption 
(A.l)' and (|4J0)) (for X) that 

E[A^„"]^A^\"[„,Z] = a^(t)E[A^-XA^\"[„,X]+0(n-^(*)-i+n-«(*)-) 

implying (I4.10p for Z (with c{t) replaced by a'^{t)c{t)). Whence and using (|4.10p and (|6.18p . relation (I6.20p 
follows similarly as (j6.16l) above. 

The proofs of (|6.2ip - (|6.22p uses the following bounds from Bruzaite and Vaiciulis (2008, Lemma 1). Let 
(Ui, U2) G be a Gaussian vector with zero mean, unit variances and a correlation coefficient p, \p\ < 1. 
Then for any 61, 62 G K and any 1/2 < ai, 02 < 2 

|EV'(aiJ7i+fei,a2C/2 + &2)-EV'(aiJ7i,a2t/2)| < C{bl+bl), (6.23) 
\EU,i:{aiUi + bi,a2U2 + b2)\ < C(|6i| + I&2I), « = 1,2, (6.24) 
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where the constant C depends only on p and does not depend on 01,02,61,62. Using (|6.23p . Assumption 
(A.l)' and the fact that |A^'"/3| < Cn"^ we obtain 

n — 3 

|Ei?"'2 _E^"'2| < C'7i-i^ (n^e^/")?!-^)^ = 0{n-^), (6.25) 

fc=0 

proving (|6.2ip . 

To prove O, write i?"'^ ^ _j_ ^"-3 ^^(^^^ ^„,2 ^ _i_ ^j^^^^ 

r;„(fc) := ^(A^-"Z,A^:;\Z) = /„,fe(yW(fc) + MW(fc),yi2)(fc)+Ml')(fc)) 

where standardized increments Y,i^\k),i = 1,2 arc defined as in ()6.9p - ()6.10p with X replaced by Z, fn^k are 
defined in (|6.11l) . and 



Note, the 77„(fc)'s and ?7„(fc)'s have the Hermite rank > 2 and > 1, respectively, since the Hermite coefficients 
of order 1 ciP {k) E[ri*^(fc)7?„(fc)] , i = 1,2 of the 77„(fc)'s are not zero in general. Using the bomrd in 
(lOi)) and |A^'"^| < Cn-^ we obtain 

|c«(fc)| < C{U'Hk)\ + \^^^'\k)\) < Cn^C^/")-^ (6.26) 

Split r7„(fc) - 77„(fc) - (E?7„(fc) - Ery„(fc)) = Caik) + C{k), where 
2 

C(fc) := ^c«(fc)rW(fc), C(fc) := Vn{k) - Cik) - Uk) ~ (E?7„(fc) - Ery„(fc)). (6.27) 

4=1 

Then var(i?"'2 _ j^n,2) < 2(j; + j;;), where J^, := e(;A_ ^;(fc)) j;; := e(;A_ ^;'(fc)) 

From (IQej) and Cauchy-Schwartz inequality, it follows j; < (;^^((n - 2) ^^I-Iq n^-f^^*^/")-^^ ^ 0(l/n2) = 
o(l/n). Since the have Hermite rank > 2, we can use the argument in the proof of Theorem 14.21 

together with Arcones' inequality (Arcones, 1994) and the easy fact that E{(^{k))'^ — ^ (n — > 00, k = 
0, • • • , n — 3), to conclude that J" = o{l/n). This proves (|6.22p and completes the proof of part (i). 

(ii) Similarly as in the proof of part (i), we shall restrict ourselves to the case p = 2 for concreteness. We 
use the same notation and the decomposition of i?^'" as in part (i): 

i?"'" - Ap{H{T))dT = - Ap(iJ(r))dT) + (Ei?"-f - ER"'P) + Q„, (6.28) 

where 

n — 3 n — 3 

see ()6.27p . Here, the a.s. convergence to zero of the first term on the r.h.s. of (|6.28p follows from Theorem 
14. II since the conditions of this theorem for Z are easily verified. The convergence to zero of the second term 
on the r.h.s. of (|6.28p follows similarly as in (|6.25p . with the difference that |A^'"/3| < Cn~^ since /3 € C^[0, 1] 
and therefore 

n—3 

|Ei?"^2 _E^«,2| < Cn-i^ (n^('=^/")n-i)^ = 0(n~2(i-sup,^[„_ij ^ ^^^^^ ^g_29) 

k=0 

The proof of Q'„' ^ mimics that of Theorem 14.11 and rehes on the fact that the have the Hermite 

n— >-oo 

rank > 2 (see above). Relation Q'^ follows by the gaussianity of Q'^ and E((5^)^ — 0{n^^) for some 

S > 0. Since |cL'^(fc)| < Cn^C^/")-! ^^g^ ^^^^ inequality in (lOe]) and |A^'"/3| < Cn^\ then 

E((5'J2 < Cn~^(^"™P*e[o.ii -^('^^ similarly as in (j6.29p above. This proves part (ii) and the proposition. □ 
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Proof of Corollary 14- 2 



(a) The argument at the end of the proof of Property 4.1 shows that V satisfies Assumption (A.l)', while 
(A.2)p follows from Property 14.21 Then the central limit theorem in (14.221) follows from Theorem 4.2. 

(b) In this case, (A.2)p follows from Property 14.31 Instead of yerifying (A.l)', it is simpler to directly yerify 
condition (|4.10|) which suffices for the validity of the statement of Theorem 4.2. Using /(^) = c^~^^~^(l + 
o(^"^/^)) — >■ oo), similarly as in the proof of Property 14.31 for j e N* one obtains 



poo 

^ n-2^^E[Ag'"XAf"X] -c2i+P / (1 - cos(x))f cos(xj) : 

Jo 



-2H-1 



dx 



2i+P 



(1 ~ cos{x)f cos{xj) X ^ {n^"+^ifinx) - c{nx)-^"-^))dx 



by the Lebesgue dominated convergence theorem, since |(1 — cos{x)y cos(xj) x^^^^"^^^\dx < oo. There- 
fore condition (|4.10p is satisfied and Theorem 14.21 can be applied. 

Finally, the function H K2{H) is a C^(0, 1) bijective function and from the Delta-method (see for instance 
Van der Vaart, 1998), the central limit theorem in (|4.23p is shown. □ 



Proof of Lemma \5.1\ 

(i) We use the following inequality (see, e.g. Ibragimov and Linnik, 1971, Theorem 1.5.2): 



\Fn — GqIIoo < 



2tt 



fn{0)-gc.iO) 



d0. 



(6.30) 



where /„, ga are the characteristic functions of Fn, Ga, respectively. According to (|5.6p and the definition of 

fn{e)-g^{e) = g„(0)(e-2n-^/.W_l)^ 

where, with v :— 9n^^°', 



l/v 



Ini0) := / Re(e'"''-^l)di^i(M) = 

Jo Jo 

If w > 1 then integrating by parts and using (|5.9p . 



l/v 



I/2I < 2|ifi|(l/«) = 0(«("-'5)+), 

\Ii\ = A'i(l/w)Re(e' - 1) - / Ki{u) v sm{uv) du 

Jo 

< 2\Ki{l/v)\+v^ u\Ki{u)\du 
Jo 

Next, if u < 1, then similarly as above 

I/2I < 2|i^i|(l/t;) = 0{v'), 



Ki{l/v)Re{e -I) - Ki{u) V siii{uv) 
Jo 



du 



l/v 



< 2\Ki{l/v)\+v' / u\Ki{u)\du 
Jo 
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Therefore, for some constant C, 

Moreover, gaiO) — e^'^^'^l for some ci > 0. Using these facts and ()6.30p we obtam 

||F„-G„||.o < ci f'"' n-\en''n''''^ + + r n-i(0ni/")("-^)+e-^^l«l°^ 

1 

The above bound easily yields (|5.10p . 

(ii) Follows similarly using ()6.30p and the argument in (i) with Ki replaced by K. □ 
Computation of X{r) 

From the definition of A(r) and the change of variables Xi = acoscf), X2 = asin0, with |r| < 1, 

1 f \X1+X2\ —^^(x\-2rx^x^^xl) 



COS0 + sm0| 

d0 



TT 7o (I COS 01 + I sin — r sin(20)) 
/1+/2, 



where 



^ n Jq l~r sin(20) 



VT^^ /""^ 1 , 1 1 / r \ 1 

-di = — ^ arc.ta.n = — arccosi— r ; 



TT Jo l+t^-2rt 2 TT V^l _ 7.2/ TT 

2\/l~72 cos0-sin(/) ,, VT"^ /"i 1-i 



dt 



lo: 



(cos(/i + sin(/))(l + rsin(20)) Jo (l + t){l + 2rt + 1^) 

7r(l - r) Vr + l/ 

The function A(r) is monotone increasing on [—1, 1]; A(l) = 1, A(— 1) = 0. It is easy to check that 

o2H I n2H+2 _ n 

p,{H) = 2^^-1-1, p2{H)^ 3_^,g+i (6.31) 

are monotone increasing functions; pi(l) = l,p2(l) = so that Ap(iJ) = X{pp{H)) for p = 1,2 is also 
monotone for e (0, 1). 

Expression and graph of S2{H) 

[ d 1 ^ 

From the Delta-method, sl{H) = — (A2)"^(A2(H)) Y^2{H) and therefore 



2,^.( ^(8- 2^^+1)2(1 -p2(g))Vl-p^(H) 

' V(log2-log(l+p2(i?)))(22«+291og2-32«161og3 + 62«41og(3/2))/ ^' 

with the approximated graph (using the numerical values of J^2{H) in Stoncelis and Vaiciulis (2008)): 
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